{
 "cells": [
  {
   "cell_type": "code",
   "execution_count": 2,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Populating the interactive namespace from numpy and matplotlib\n"
     ]
    }
   ],
   "source": [
    "# referecence: NCBI website\n",
    "%pylab inline\n",
    "from bioservices import EUtils\n",
    "from easydev import AttrDict"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 4,
   "metadata": {},
   "outputs": [],
   "source": [
    "e = EUtils()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 5,
   "metadata": {
    "collapsed": true
   },
   "outputs": [],
   "source": [
    "report = e.ESearch('nuccore', term='A/chicken/Hong+Kong/915/97+OR+A/chicken/Hong+Kong/915/1997', usehistory='y')"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 6,
   "metadata": {
    "collapsed": true
   },
   "outputs": [],
   "source": [
    "report = AttrDict(**report)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 8,
   "metadata": {},
   "outputs": [],
   "source": [
    "import pandas as pd\n",
    "summary = e.ESummary('nuccore', WebEnv=report.webenv, query_key=report.querykey)\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 9,
   "metadata": {
    "collapsed": true
   },
   "outputs": [],
   "source": [
    "uids = summary['uids']"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 10,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/html": [
       "<div>\n",
       "<style>\n",
       "    .dataframe thead tr:only-child th {\n",
       "        text-align: right;\n",
       "    }\n",
       "\n",
       "    .dataframe thead th {\n",
       "        text-align: left;\n",
       "    }\n",
       "\n",
       "    .dataframe tbody tr th {\n",
       "        vertical-align: top;\n",
       "    }\n",
       "</style>\n",
       "<table border=\"1\" class=\"dataframe\">\n",
       "  <thead>\n",
       "    <tr style=\"text-align: right;\">\n",
       "      <th></th>\n",
       "      <th>accession</th>\n",
       "      <th>description</th>\n",
       "      <th>ids</th>\n",
       "      <th>taxid</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>0</th>\n",
       "      <td>AF098594.1</td>\n",
       "      <td>Influenza A virus (A/chicken/Hong Kong/915/97(...</td>\n",
       "      <td>gi|6048875|gb|AF098594.1|</td>\n",
       "      <td>100838</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1</th>\n",
       "      <td>AF098581.1</td>\n",
       "      <td>Influenza A virus (A/chicken/Hong Kong/915/97(...</td>\n",
       "      <td>gi|6048849|gb|AF098581.1|</td>\n",
       "      <td>100838</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2</th>\n",
       "      <td>AF098550.1</td>\n",
       "      <td>Influenza A virus (A/chicken/Hong Kong/915/97(...</td>\n",
       "      <td>gi|6048770|gb|AF098550.1|</td>\n",
       "      <td>100838</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>3</th>\n",
       "      <td>AF098564.1</td>\n",
       "      <td>Influenza A virus (A/chicken/Hong Kong/915/97(...</td>\n",
       "      <td>gi|6048802|gb|AF098564.1|</td>\n",
       "      <td>100838</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>4</th>\n",
       "      <td>AF098620.1</td>\n",
       "      <td>Influenza A virus (A/chicken/Hong Kong/915/97(...</td>\n",
       "      <td>gi|6048927|gb|AF098620.1|</td>\n",
       "      <td>100838</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>5</th>\n",
       "      <td>AF098608.1</td>\n",
       "      <td>Influenza A virus (A/chicken/Hong Kong/915/97(...</td>\n",
       "      <td>gi|6048903|gb|AF098608.1|</td>\n",
       "      <td>100838</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>6</th>\n",
       "      <td>AF098573.1</td>\n",
       "      <td>Influenza A virus (A/chicken/Hong Kong/915/97(...</td>\n",
       "      <td>gi|6048829|gb|AF098573.1|</td>\n",
       "      <td>100838</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>7</th>\n",
       "      <td>AF046100.1</td>\n",
       "      <td>Influenza A virus (A/chicken/Hong Kong/915/199...</td>\n",
       "      <td>gi|3421265|gb|AF046100.1|</td>\n",
       "      <td>380304</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "</div>"
      ],
      "text/plain": [
       "    accession                                        description  \\\n",
       "0  AF098594.1  Influenza A virus (A/chicken/Hong Kong/915/97(...   \n",
       "1  AF098581.1  Influenza A virus (A/chicken/Hong Kong/915/97(...   \n",
       "2  AF098550.1  Influenza A virus (A/chicken/Hong Kong/915/97(...   \n",
       "3  AF098564.1  Influenza A virus (A/chicken/Hong Kong/915/97(...   \n",
       "4  AF098620.1  Influenza A virus (A/chicken/Hong Kong/915/97(...   \n",
       "5  AF098608.1  Influenza A virus (A/chicken/Hong Kong/915/97(...   \n",
       "6  AF098573.1  Influenza A virus (A/chicken/Hong Kong/915/97(...   \n",
       "7  AF046100.1  Influenza A virus (A/chicken/Hong Kong/915/199...   \n",
       "\n",
       "                         ids   taxid  \n",
       "0  gi|6048875|gb|AF098594.1|  100838  \n",
       "1  gi|6048849|gb|AF098581.1|  100838  \n",
       "2  gi|6048770|gb|AF098550.1|  100838  \n",
       "3  gi|6048802|gb|AF098564.1|  100838  \n",
       "4  gi|6048927|gb|AF098620.1|  100838  \n",
       "5  gi|6048903|gb|AF098608.1|  100838  \n",
       "6  gi|6048829|gb|AF098573.1|  100838  \n",
       "7  gi|3421265|gb|AF046100.1|  380304  "
      ]
     },
     "execution_count": 10,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "# here you need a little bit of introspection but can easily obtain the accessions, descriptions, ids\n",
    "accession = [summary[s]['accessionversion'] for s in uids]\n",
    "description = [summary[s]['organism'] for s in uids]\n",
    "id_ = [summary[s]['extra'] for s in uids]\n",
    "taxid = [summary[s]['taxid'] for s in uids]\n",
    "\n",
    "\n",
    "pd.DataFrame({'accession':accession, 'description': description, 'ids':id_,'taxid':taxid})\n",
    "    "
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 11,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "b'>AF098594.1 Influenza A virus (A/Chicken/Hong Kong/915/97 (H5N1)) PB1 protein (PB1) gene, partial cds\\nGTCAATCCGACTTTACTTTTCTTAAAAGTGCCAGCGCAAAATGCCATAAGTACCACATTCCCTTATACTG\\nGAGATCCCCCATACAGCCATGGAACAGGAACAGGGTACACCATGGACACAGTCAACAGAACACATCAATA\\nTTCAGAAAAAGGGAGGTGGACAACAAACACAGAGACCGGAGCACCCCAACTCAACCCTATTGATGGACCA\\nTTACCTGAAGATAACGAGCCGAGCGGGTATGCACAGACAGATTGTGTATTGGAAGCAATGGCTTTCCTTG\\nAAGAATCCCACCCAGGACTCTTTGAAAACTCATGTCTTGAAACGATGGAAGTTGTCCAGCAAACGAGAGT\\nGGATAAGCTGACCCAAGGTCGCCAGACCTATGACTGGACATTGAACAGAAACCAGCCGGCTGCAACTGCT\\nTTGGCCAACACCATAGAAGTATTCAGATCGAACGGTCTAACAGCCAATGAGTCAGGAAGGCTAATAGATT\\nTCCTCAAGGACGTAGTGGAATCAATGGATAAGGAAGAAATGGAGATAACAACACATTTCCAGAGAAAGAG\\nAAGAGTGAGGGACAACATGACCAAGAGAATGGTCACACAAAGAACAATAGGGAAAAAGAAGCAAAGGCTG\\nACAAAAAAGAGCTACCTAATAAGAGCACTGACACTGAACACAATGACAAAAGATGCTGAAAGGGGAAAAT\\nTGAAAAGGCGAGCAATTGCAACACCCGGAATGCAAATCAGAGGATTCGTGCACTTTGTCGAAGCACTAGC\\nAAGGAGCATCTGTGAAAAACTTGAGCAATCTGGACTCCCCGTCGGAGGGAATGAGAAGAAGGCTAAATTG\\nGCAAATGTTGTGAGAAAGATGATGACTAACTCACAAGATACAGAGCTCTCCTTTACAGTTACCGGAGACA\\nACACCAAATGGAATGAGAATCAGAATCCTCGAATGTTTCTAGCAATGATAACATACATCACAAGGAACCA\\nACCTGAATGGTTTAGAAATGTCTTGAGCATTGCCCCTATAATGTTCTCAAATAAGATGGCGAGGTTAGGG\\nAAAGGATACATGTTCGAGAGTAAGAGCATGAAGCTACGGACACAAATACCAGCAGAAATGCTTGCAAACA\\nTTGACTTGAAATACTTCAACGAATCAACGAGAAAGAAAATTGAGAAAATAAGACCTCTACTAGTAGAGGG\\nCACAGCCTCATTGAGTCCAGGGATGATGATGGGCATGTTTAATATGCTAAGTACGGTCTTAGGAGTCTCA\\nATCTTAAATCTTGGGCAGAAGAGGTACACCAAAACCACATACTGGTGGGATGGGCTCCAATCCTCTGATG\\nATTTCGCTCTCATAGTGAATGCACCAAATCATGAGGGAATACAAGCAGGAGTGGATAGATTCTATAGGAC\\nTTGCAAACTAGTTGGAATCAACATGAGCAAAAAGAAGTCTTACATAAATCGGACAGGAACATTTGAGTTC\\nACAAGCTTTTTCTACCGCTATGGGTTTGTAGCCAACTTCAGCATGGAGCTGCCCAGCTTTGGAGTTTCCG\\nGAATTAATGAATCGGCTGACATGAGCATTGGAGTTACAGTGATAAAGAATAATATGATAAACAACGACCT\\nTGGACCAGCAACAGCTCAGATGGCTCTTCAGCTGTTCATTAAGGACTACAGATACACCTATCGATGCCAC\\nAGAGGTGATACACAAATTCAAACTAGAAGATCATTTGAATTGAAGAAGCTGTGGGAGCAGACCCGCTCAA\\nAGGCAGGACTGTTGGTTTCAGATGGAGGGCCGAATTTATACAACATCCGGAATCTTCACATTCCAGAAGT\\nTTGCTTGAAGTGGGAGTTGATGGATGAAGATTATCAGGGAAGACTGTGTAACCCTCTGAACCCGTTTGTC\\nAGTCATAAGGAAGTTGAATCCGTCAACAATGCTGTGGTAATGCCAGCCCATGGTCCAGCCAAGAGCATGG\\nAATATGATGCCGTTGCAACTACACATTCATGGATTCCCAAGAGAAACCGCTCCATTCTCAACACTAGCCA\\nAAGGGGAATTCTTGAGGATGAACAAATGTACCAGAAGTGCTGCACTCTATTCGAGAAATTCTTCCCTAGC\\nAGTTCATACCGGAGGCCAGTTGGAATTTCCAGCATGATGGAGGCCATGGTGTCTAGGGCCCGAATTGATG\\nCACGGATTGACTTCGAATCTGGAAGGATTAAGAAAGAAGAATTTGCTGAGATCTTGAA\\n\\n'"
      ]
     },
     "execution_count": 11,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "sequence = e.EFetch('nuccore', uids[0], rettype='fasta', retmode='text', WebEnv=report.webenv, query_key=report.querykey)\n",
    "sequence"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "### Using EGQuery to get H5N1 Related-Records in Entrez\n",
    "\n",
    "To see what other Entrez databases contain information about the H5N1 virus, use EGQuery. EGQuery performs a text search across all available Entrez databases and returns the number of hits in each. EGQuery accepts any valid Entrez text query as input through the term parameter."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 12,
   "metadata": {
    "collapsed": true
   },
   "outputs": [],
   "source": [
    "results = e.EGQuery(\"H5N1 virus\")"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 13,
   "metadata": {
    "collapsed": true
   },
   "outputs": [],
   "source": [
    "\n",
    "data = [x['Count'] for x in results['eGQueryResult']['ResultItem'] if int(x['Count'])>0]\n",
    "index = [x['DbName'] for x in results['eGQueryResult']['ResultItem'] if int(x.Count)>0]\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 14,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAbAAAAEICAYAAAA+16EyAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAIABJREFUeJzt3Xm8VVX5x/HPlyFQwCnIcEDMcEIFFTAUx8g0LckhckpM\n5We/HPuhaaWhZmpa5pADmqI5K5JTKqYCzkwyOQ9glmSOICqp8Pz+WM+BzeGcc8+93HvOPfc+79fr\nvu65e6+99jr7DuuuvdezHpkZIYQQQq1pU+0GhBBCCA0RHVgIIYSaFB1YCCGEmhQdWAghhJoUHVgI\nIYSaFB1YCCGEmhQdWKgoSaMl/aZK55akayV9IGlSNdrQHEm6X9Jhzb3OlSFppKQbVuL48ZKObMw2\nZequ2u9ErYsOrJWTNFfS25I6ZbYdKWl8FZvVVAYB3wLWM7MB+TslDZP0eIHtcyUN9te7SFoiaWHm\n47C8siWvp6SzJM2S9IWkkY37Fpdrd1l/tM1sTzO7rjHP3RR11orsz0toWtGBBYB2wPHVbkR9SWpb\nz0M2AOaa2ccreeq3zKxz5iP/D3Vd1/NV4GTgvpVsx0rxEWn8DQg1K354A8D5wAhJa+TvkNRTkklq\nl9m29HaKj1qekHShpA8lvS5pe9/+pqT/FLiV1FXSQ5I+kjRB0gaZujf1fe9LeknSDzL7Rku6XNLf\nJH0M7FqgvetIutuPf1XSUb79COBqYKCPms5YyWtWStHrCWBm15nZ/cBHdVUkqY2kUyS9Juk9SbdJ\nWsv35b43h0n6h6R3Jf3S9+0B/AIY6u93hm8fL+lsSU8AnwBfy789JunHkl7wW60P5r4/3uFd6N/T\n+ZJmStqiSLvzf0Yel3SB1zlH0p4l3vM6ksZIesfLHpfZN0DSU/6zNk/SpZK+lNnfO/Pz87akX2Sq\n/pKk6/3n7jlJ/Uq04VuSXvT3eSmgzL6NJD3i3493Jd2Y+15L+gvQA7jHr/vJvv12Sf/2+iZK6p13\nylK/Exf579ICSVMl7Zh3Pab4vrcl/SGz7xuSnvRrNUPSLsXeb62KDiwATAHGAyMaePx2wEzgy8BN\nwC1Af+DrwCHApZI6Z8ofDJwFdAWmAzcCKN12e8jr+ApwIHBZ3i/7QcDZQBdghdt9wM3AP4F1gP2B\n30r6ppn9GTgaeMpHTb9u4HsF+Ir/sZjjf9A75e1f2euZdRwwBNiZ9J4+AP6UV2YQsAnwTeB0SZuZ\n2QPAb4Fb/f32yZQ/FBhOuoZvZCuSNITU8e0LdAMeI11TgN2BnYCNgTWAocB7Zb6P7YCXSN/z3wF/\nlqT8QkojwnuAGcC6/p5OkPRtL7IYONHrGej7/9eP7QL8HXiAdK2+Djycqf57pJ/NNYC7gUsLNVRS\nV2AM8Cs/z2vADtkiwDl+js2A9YGRAGZ2KPAP4Lt+3X/nx9wP9CL9XE/Df+YzCv5OuMlAX2At0u/G\n7ZI6+r6LgIvMbDVgI+A2fw/rkkb4v/HjRgBjJHUr9J5rlpnFRyv+AOYCg4EtgPmkP1pHAuN9f0/A\ngHaZY8YDR/rrYcArmX1bevm1M9veA/r669HALZl9nUl/lNYn/UF8LK99VwK/zhx7fYn3sr7X1SWz\n7RxgdKatj5c4fhjwBfBh3scSYLCX+SqwOemfvw2BicCV5V7PvPPdAIys4/vzAvDNzNfdgc9Jtylz\n35v1MvsnAT/01yOBG/LqGw+cWWBb7vt5P3BEZl8b0khtA2A34GXgG0CbOtqd/zPyambfqt7urxY4\nbjvgH3nbTgWuLXKeE4Cx/vpA4Nki5UYCf898vTnwaZGyPwKeznwt0j9FRxYpPyR73tzPQIlrs4a/\n/9Xr+p0ocvwHQB9/PRE4A+iaV+bnwF/ytj0IHFbq+1ZrHzECCwCY2WzgXuCUBhz+dub1p15f/rbs\nCOzNzHkXAu+T/pvdANjOb3l8KOlD0n+mXy10bAHrAO+bWfbW3Buk/+TL9bSZrZH9IP1HnWvvv83s\neTNbYmZzSM+y9s+vZCWvZ9YGwNjM9XiB9Mdt7UyZf2def8Ly17qQUtdwA+CizPneJ/0BX9fMHiGN\nWv4EvC1plKTVynwfS9toZp/4y0Lt3ABYJ+9n4Bf4+5W0saR7/XbcAtIos6sfuz5ptFRnG0jXqaMy\nt8Yz1mH5n1HLfi3pK5JukfQvb8MNmTasQFJbSecq3QZeQOrgyDum2O8Ekv7Pb+nO9+uxeubYI0gj\n4hclTZa0t2/fADgg7zoOIv0D1GJEBxayfg0cxfJ/8HMTHlbNbMt2KA2xfu6F31pcC3iL9Es8Ia8D\n6WxmP8kcWyp9wlvAWn4rKacH8K+VbG8pRub5SJ5C17O+3gT2zLsmHc2snPdU7FqVuoZvAv+Td75V\nzOxJADO72My2BXqT/nCeVJ83U4Y3gTl55+9iZt/x/ZcDLwK9LN02+wXLrv+bpNtoK2sey/+MKvs1\naVRvwFbehkNY/mcg//oeBOxDGpmvTho5k3dMwd8Jf971c+AHwJr+D9X83LFm9oqZHUi6NXkecIff\n0n6TNALLXsdOZnZufS9GcxYdWFjKzF4FbiU9d8lte4fUARzi/0n+mJX/I/EdSYP84ftZwDNm9iZp\nxLKxpEMltfeP/pI2K7P9bwJPAudI6ihpK9J/qPnPGxpMaRp9DyXrA+cCdxVpzwrX0+to788w2gDt\nvK3FZlReAZytZRMpuknap8zmvg30VP1mGl4BnJp77ihpdUkH+Ov+kraT1J70j80i0miwMU0CFkj6\nuaRV/GduC0n9fX8XYAGwUNKmQPafm3uBr0o6QVIHSV0kbdeANtwH9Ja0r4/QjmP5f9q6AAuBD/1Z\nU34n/jbwtbzy/yXdSl+VNGrMV+x3ogvptvY7pJ+V04Glo15Jh0jqZmZLSLe7IX1PbgC+K+nbfg07\n+s/uevW/HM1XdGAh35lA/qSEo0i/pO+R/vN+ciXPcRNpdPI+sC3pNiF+62934Iek0dS/Sf9VdqhH\n3QeS/sN9CxhLen720Eq2N2sb4CnSH/AngdnkdVB5Cl3Pq0i3VQ8EfumvDy1y/EWkCQfjJH0EPE16\nTlSO2/3ze5KmlXOAmY0lXfNb/HbXbCA3Y3A1b/sHpFuz7wEXlNmWspjZYuC7pEkLc4B3SbNHV/ci\nI0gjmo+8Lbdmjv2IFOf3XdLPzisUmKlaRhveBQ4g/XPyHmnyxROZImeQfg7mkzq7O/OqOAf4ld+6\nGwFcT7pe/wKeJ30P8xX8nSA9t7qf9OzxDdI/DdlbwHsAz0laSPpZ+aGZLfLObx/SCPUdP+YkWtjf\nfPnDvRBCCKGmtKjeOIQQQusRHVgIIYSaFB1YCCGEmhQdWAghhJpUKIgvNJKuXbtaz549q92MEEKo\nKVOnTn3XzOpc9qrVdWCSzge+A/yNFLX/iZldL2k0cK+Z3dFY5+rZsydTpkxprOpCCKFVkPRG3aVa\nYAcmqZ2ZfVGiyP8A3czsv03dlln/mk/PU6qaMSOEECpu7rl7VeQ8FXsGJuk0pfQED0m6WdIIpbQE\nDyilCHjMI+uRtIGkh5XSNTwsqYdvHy3pCi/7cm7dL6V0DbdLuocU8ClJ50uarZQ4cKiXu5sUVPqM\npKFKCf9WWDFc0rZKKQ2mKqWT6O7bj5P0vLfrlgpduhBCCAVUZASmlHdnP2BrP+c0YCowCjjazF7x\nJV8uI614fSlp1fHrfOmii0krPkNaZWFn0nJGj0r6um8fSFqb7H1J+5Ei+fuQFr2cLGmimX1P0kIz\n6+vtGlmgre2BS4B9zOwd7/zOBn5MWph1QzP7r4rkepI0nJSqgrartazMBSGE0JxU6hbiIOAuM/sU\nwEdKHYHtSbltcuVySwYNJOUjAvgLKX9Qzm2+7tcrkl4HNvXtD5nZ+5nz3ezL0rwtaQIpP9XdZbR1\nE1IqjIe8XW1Ji3tCynl1o6S/An8tdLCZjSJ1zHTo3iuWOQkhhCZSqQ6s0GrdbYAPc6OhOliR19mv\ns2nii60OXg4Bz5nZwAL79iIl9PsecJqk3qWet2257upMqdC94BBCaG0q9QzscdLKyB2VUgXsRcrH\nMyez0rUk5bLGPkla0BXSopbZzLsHKKVZ34i04vNLBc43kZRKva1SBtKdSKtcl+MloJukgd6u9kpp\nytuQEsw9SsoBtQZ1510KIYTQRCoyAjOzyT6BYgZpReUppJWcDwYul/QroD0p3fcM0ure10g6ibSS\n8uGZ6l4CJpAS3B1tZou0YmbysaTbkDNII7STzezf+YWKtPUzSfsDF0tanXSN/khaDfoG3ybgQjP7\nsERVIYQQmlDFVqOX1NnMFkpalTRCGm5mZaV4yNQxmkaI1ZI0DBhnZm/VUe5MYKKZ/b0h5+nXr59F\nHFgIIdSPpKlm1q+ucpWMAxslaXPS5I3r6tt51Zektj6Jo5BhpDxHJTswMzt9ZdoQcWAhhFpUqTiu\nlVWxDszMDmqEOoYBSOoJPAA8Q5qa/zLwI1KyuGtISREvlfQiKcPsqqRVN34MfBPoR5pN+CnpVuPm\nwB9Iz7TeBYaZ2bzsiE/SXOA6UrK89sABZvbiyr6nEEIIDVPLi/luAowys61IKcb/17cvMrNBZnYL\nKRPqz73MLFJ23jtIz+AO9hmQX5DivvY3s21JHeDZRc75rpltA1xOygy7AknDJU2RNGXxJ/Mb552G\nEEJYQS0vJfWmmeXSfN/AsrTutwL4ZIs1zGyCb7+OZSnWs0rFfeXLpQ6fyrI4teVEHFgIIVRGLXdg\n5cSDlaNU3Fe+3PqJiynj2kUcWAghNJ1avoXYIxerBRzI8rFimNl84ANJO/qmQ0nT7wE+Arr464Jx\nX03a8hBCCCutljuwF4DDJM0E1iI9l8p3GHC+l+kLnOnbRwNXSJpOumW4P3CepBnAdNISVyGEEJqx\nisWBNSafhXivmW1R5aaUFHFgIYRQf80xDqxmlZFjrKCIAwshNFe1EutVSk3cQpTUU9ILkq6S9Bxp\nll9/SeM9VQuSunqsFr4G4gWeC2ympGN9e39JT0qaIWmSpC6+PuO1XvZZSbt62eVyjPm2kyRN9jrP\nqMa1CCGEkNTSCKwXcKCZHSXpNlJ+sWKGAxsCW5vZF5LWkvQl0hT7ob4242rAp8DxAGa2pVJCzXGS\nNvZ6sjnGdvc2DCDNXLxb0k5mNjF74sgHFkIIlVETIzA3x8ym++uppMSWxQwGrsjd9vM8YZsA88xs\nsm9b4PsHkXKO4StrvAHkOrBsjrHd/eNZUkLOTUkd2nLMbJSZ9TOzfm1XXb2h7zWEEEIdamkE9t/M\n68XAKqRVNHKdcMfMfrFinFihbbntxeTnGDvHzK4sq7UhhBCaVC11YIXMBbYl5fraP7N9HHC0pPG5\nW4jAi8A6kvr7LcQupFuIE0lpXR7xW4c9SLFh2+Sd60HgLEk3+qr66wKfm9l/ijUuAplDCKHp1NIt\nxEIuAH4i6Umga2b71cA/gJke23WQmX0GDAUu8W0PkUZtlwFtJc0iPSMbZmbZ0R4AZjYOuAl4ysve\nwbJg6BBCCBVWk3FgtSLiwEIIof7KjQOr9RFYCCGEVqrWn4HVSx15xG4CdiXl+hoOnAN8HTjfzK7w\n408mram4BLjfzE4pdb4IZA4h1KUlBBRXS6vqwNwmwBFm9oSka1iWR+xNMxso6ULSWok7kJ6RPUda\nN3FPYAiwnZl94hNDQgghVElrvIWYn0dskL++2z/PAp4xs4/M7B1gkaQ1SLFl15rZJ7A0tmwFkdAy\nhBAqozV2YMXyiOVmHi5h+ZizJaSRarE4suUri0DmEEKoiNZ4C7GHpIFm9hTL8ohtXcZx44DTJd2U\nu4VYbBSWE3FgIYTQdFrjCKycPGIrMLMHSLcZp3gesRFN18QQQgh1aXFxYJJGAgvN7IIC+3rSSHnE\nfOX7fmb2brEyEQcWQgj1F3FgIYQQWrRm/wysjtitfmb2rucEu8DMdvHD+kh6BFgf+J2ZXSVpF+AM\n4HlJLwN3kmYcHk9aGHiImb0mqRtwBWlNRIATfMr9l4GbgW6ktRdLLQIMRBxYCK1FxHJVR62MwDYB\nRpnZVsAClsVuFbMVsBcpn9fpktbx7X1IHdaWpIDkjc1sAGntxGO9zEXAhWbWn5Rz7Grf/mvgcTPb\nmvQsLNfBhRBCqIJmPwJz+bFbx9VR/i4z+xT4VNKjpCSUHwKTzWwegKTX8EzLpJHYrv56MLC5tHSA\ntZqvXL8TsC+Amd0n6YNCJ46EliGEUBm10oEVit0qlgusWHlYMb4rG/uVuxZtgIHeAS7lHVpZcWDA\nKIAO3Xu1rBkyIYTQjNRKB1YodqsLKRfY/aRbfVn7SDoH6ATsApzCsizLdRkHHAOcDyCpr2eCzuUN\n+40vK7VmXRVFHFgIITSdWnkGVih26wzgIkmPkTI0Z00C7gOeBs4ys7fqca7jgH6SZkp6Hjjat58B\n7CRpGrA7Kd9YCCGEKmn2cWCNGbtVaREHFkII9RdxYCGEEFq0Zv8MzMzmAk0y+pJ0Gum51pvAu8BU\nYD5pFuGXgFeBQ33tw9HAIqA3sDbwMzO7t1T9EQcWQt0ihio0VKsdgXnw836k4Oh9gdxw9U4z629m\nfUjP3o7IHNYT2JkUY3aFpPzZjyGEECqk1XZgpDxgd5nZp2b2EXCPb99C0mOSZpFGZ70zx9xmZkvM\n7BXgdWDT/EojH1gIIVRGa+7Aii0FNRo4xsy2JM08zI6yisWXLdsQ+cBCCKEimv0zsCb0OHClx4u1\nI90WvIoUXzZPUnvSCOxfmWMOkHQdsCHwNeClUieIOLAQQmg6rbYDM7PJku4GZgBvAFNIEzhOIy0c\n/AZpiakumcNeAiaQJnEcbWaLKtroEEIIS7XaDsxdYGYjJa1KWmnj92Y2jeJJLp8wsxMr17wQQgjF\ntPYObJSkzUnPua7zziuEEEINaPYrcdSyDt17WffD/ljtZoRQNRHjFRqixazEIamnpBclXefrE94h\naVVJ/SU9KWmGpEmSunjZxyRN84/tvY7ukiZKmi5ptqQdJbWVNNq/niXpRC97lKTJXu8YP1cXSXN8\nYgeSVpM0N/d1CCGEymv2HZjLT2h5DHArcLwHHA8GPgX+A3zLzLYBhgIX+/EHAQ+aWV9SUsvpQF9g\nXTPbwqfMX+tlVwhk9jix8aSZigA/BMaY2ef5DY04sBBCqIxa6cDyE1p+G5hnZpMBzGyBmX0BtAeu\n8iDk24HN/ZjJwOGSRgJbeof0OvA1SZdI2oPUMULxQOargcP99eEs6/CWE3FgIYRQGbUyiSP/Qd0C\noEOBcicCb5NGWW1IaxdiZhMl7UQaQf1F0vlmdr2kPqTO8KfAD4AfkwKZh5jZDEnDSPnEMLMn/Bbl\nzkBbM5tdV6MjDiyEEJpOrYzAekga6K8PJOX5WkdSfwB/RtUOWJ00MlsCHAq09f0bAP8xs6uAPwPb\nSOoKtDGzMaTYr228/vxA5qzrgZspMvoKIYRQObUyAssltLwSeAW4BHgEuETSKqTnX4OBy4Axkg4A\nHgU+9uN3AU6S9DmwEPgRsC5wraRcJ36qfy4VyHwj8BtSJxZCCKGKmv00+vomtPRV5n9kZsc1QVv2\nB/Yxs0PLKR8JLUMIof7KnUZfKyOwspnZFNKyUI1K0iXAnsB3yj0m8oGFlixivEK11cozsKUkjZA0\nUtJ4Sed5DNjLknb0/btIutdfd5P0kMeEXSnpDX/2haRD/Njpvi/3vGyhpLM9DuxpSWv7qc8krZt4\no8eJ7VCFtx9CCMHVXAeWp52ZDQBOAH5dYP+vgUc8Lmws0ANA0makOLEdPDZsMcsmbHQCnvY4sInA\nUb79IuBCM+tPSoR5daEGRRxYCCFURq3fQrzTP08lZUvONwj4PoCZPSDpA9/+TWBbYLIkgFVIQdAA\nnwH3Zur9lr8eDGzu5QFWk9TFY8qWMrNRwChIS0k19I2FEEIorRY6sC9YfqSYTTD5X/+8mMLvpVjS\nSpEW7z21wL7PbdnMlmy9bYCBZvZpWa0m4sBCCKEp1cItxLeBr0j6sqQOwN71OPZxUoAyknYH1vTt\nDwP7S/qK71vLY8VKGUdawgo/pm892hFCCKGRNfsOzNcbPJMUm3Uv8GI9Dj8D2F3SNNIMwnnAR2b2\nPPArYJykmcBDQPc66joO6OcLCj8PHF2/dxJCCKExNfs4sJXhI7bFZvaFr+RxuU/aqIiIAwshhPpr\ntXFgeXoAt/lqG5+xbEZhvUhqa2aL63tcxIGFWhJxXaHWtOgOzMxeAbauq5ykTsBtwHqk9RPPAs4D\nrgF2By6V1AUYDnwJeBU41Mw+aaKmhxBCqEOzfwZWIXsAb5lZH1+y6gHfvsjMBpnZLRTIE1aooogD\nCyGEyogOLJkFDPaVPXY0s1zPc2umTLE8YcuJfGAhhFAZLfoWYrnM7GVJ25LWOTxH0jjf9XGm2GgK\n5AkrJeLAQgih6UQHBkhaB3jfzG6QtBAYVqBYfp6wf1WwiSGEEPJEB5ZsCZwvaQnwOfAT4I68MqXy\nhIUQQqiwFteB1Td/GICZPQg8mFdPfpnLgctXvoUhhBAaQ4vrwJqTiAMLzVXEfIWWoKXOQmwn6Tpf\n9ukOSatK+qakZyXNknSNr9JBse05klaR9ICkoyR1knSf5wqbLWlodd5eCCGEltqBbQKMMrOtgAXA\nz0izCIea2ZakkedPJHUstD1TT2fgHuAmM7uK4vFiS0UcWAghVEZL7cDeNLMn/PUNpPxfc8zsZd92\nHbATqaMrtD3nLuBaM7vevy4WL7ZUxIGFEEJltNRnYOWuUFwsX1jOE8Cekm6yZIV4MTM7s9jBEQcW\nQghNp6WOwHr46vMABwJ/B3pK+rpvOxSYQErNUmh7zunAe8BlsDRe7BMzuwG4ANimSd9FCCGEolpq\nB/YCcJjn+loLuBA4HLjdl4JaAlxhZosKbc+r6wSgo6TfkeLFJkmaDvwS+E1F3k0IIYQVtOh8YE3B\n48y2N7Ob6iob+cBCCKH+Ih9Y0+kJHATU2YFFHFiolojzCq1BTd9ClNRT0ov1iPma6zMIJ/nH1337\nAR7XNUPSRN/WVtL5kiZ73f/jpz0X2FHSdEknVuedhxBCqOkOzJUV85Upv8DMBgCXAn/0bacD3/Zc\nX9/zbUcA882sP9AfOErShsApwGNm1tfMLsxvTMSBhRBCZbSEDqzcmK+cmzOfczMVnwBGSzqKlJEZ\nUibmH/mEjWeALwO96mpMxIGFEEJltIRnYPWdhWL5r83saEnbAXsB0yX1JcWIHesL/S4laZdyTxRx\nYCGE0HRawgis3JivnKGZz08BSNrIzJ4xs9OBd4H1SavT/8TzfyFpY0mdgI+IVCohhFB1LWEElov5\nuhJ4BTgeeJoU29UOmMzysV0dJD1D6rwP9G3nS+pFGnU9DMwAZpJmHE5Tyq3yDjDEt38haQYwutBz\nsBBCCE2v2Xdgkk4gTdL4pEiRJWZ2dN62h4Gti5Rvb2bbZTeY2b4FyhnwC//I980STQ4hhFABzT6Q\nWdJcoJ+ZvVtg39eAu8tNXul1dTWzzvVsQzsz+6I+xwB06N7Luh/2x7oLhlCHiOsKrUm5gczN6hlY\ngXxbvwbWAR6V9KiXWSjpTL8NuA7QWVJX39dP0nh/3VnStR4LNlPSfsAtwCoew3Wjx5HNzpx/hKSR\n/nq8pN9KmgAcL6mbpDEeFzZZ0g4VvDQhhBDyNLdbiLl8W3sBSFqdtFbhrpkRWCdgtk+4ID2eKug0\nUhzXll5uTTMbI+kYM+vr23rW0Z41zGxnL3sTcKGZPS6pB2mSx2b5B0gaDgwHaLtat7LedAghhPpr\nbh3YLOACSecB95rZYwU6qMXAmDLqGgz8MPeFmX3QgPbcmlff5pn2rCapi5l9lD3AzEYBoyDdQmzA\nOUMIIZShWXVghfJtFSi2yMwWZ77+gmW3Qjtmtou6Y8Syx+YfD/Bx5nUbYKCZfVpHnUtFHFgIITSd\n5vYMrFC+rbriruYC2/rr/TLbxwHHZOpe019+novtAt4GviLpy75e4t4lzpNfX98631AIIYQm06w6\nMArn2xoF3J+bxFHAGcBFkh4j3V7M+Q2wZm6RXmBX3z4KmCnpRjP7HDiTtFTUvaQEl8UcB/TzCSHP\nA/lT90MIIVRQs59GD2lGIDDCzOpMrlWffF1NLfKBhRBC/bXmfGA9KTNfV1OLfGBhZUX8VwjFNdot\nxBK5ufpLetJjuyZJ6iKpt7+e7mV7lYrJcod4PbMlDfAyO3sd05Xyf3UhL1+X1/uYpGn+sb0fu4vH\net3h7b7Rl4yiSJuL5QcLIYRQBY09AtsEOMLMnpB0DWnSw9Gk3FyTJa0GfOrbLjKzGyV9iZTCZO06\n6u5kZttL2gm4BtgCGAH81M/XGVhEytc1wsz2BpC0KvAtM1uktN7hzUBuaLo10Bt4i5RSZQdJk0jT\n5/PbvDQ/mE/4eELSODObk21kxIGFEEJlNPYkjvzcXN8G5pnZZAAzW+BLMj0F/ELSz4ENypyafrPX\nMZEUg7UGqdP5g6TjSEHHhZZ7ag9cJWkWcDuweWbfJDP7p5ktAaaTbj9uUqTNZeUHi3xgIYRQGY09\nAsufEbIA6LBCIbObfCmovYAHJR0JvEzpmKz8us3MzpV0Hylu7GlJgwu06UTSdPk+Xv+izL7/Zl4v\nJl2PYvFjBfODlRJxYCGE0HQaewSWn5vraWAdSf0B/FlSO6VFeF83s4uBu4GtqDsma6jXMYh0K2++\nUh6vWWZ2HjAF2JQV48ZWJ42olpByg7WltBcLtZni+cFCCCFUQWOPwPJzc10CPAJcImkV0rOkwaTO\n6BBJnwP/Bs40s88l5WKy5rBiTNYHkp4EVgN+7NtOkLQrafT0PHA/sIRMvi7gMmCMpAOAR1l+dY0V\nmNlnkoYWaPPVFM4PFkIIoQoaLQ7M46/uLTe1SS1Y2fcUcWAhhFB/rTkOrNmIOLBQjoj1CqFhGq0D\nM7O5pKntNUPSacDBwJvAu8BU0m3Ga4BPgMczZXsD1wJfIj073M/MXql0m0MIISTNbS3EipHUj7T4\n79bAviyLDbsWOM7MBuYdkotd6+tl/1mk3uGSpkiasviT+U3T+BBCCK23AwMGAXeZ2aee0+seUrLM\nNcxsgpde/yHoAAAbo0lEQVT5S6Z8WbFrEQcWQgiV0ZqfgRVK5fwxRXKIFYpdM7NHSp0g4sBCCKHp\ntOYR2OPAdyV19GWocj3NfI81g/R8DIAisWshhBCqpNWOwHydw7uBGcAbpEDo+cDhwDWSPiEFL+es\nELtW4SaHEELIqMoILH/l+cz2qyVtXuiYJnKBmW1CCkjeBJhqZlPNrI9P4ugJjAQws3PMrLeZ9TWz\nPczs/Qq2M4QQQp5mNQIzsyMrfMpR3mF2BK4zs2mNWXnEgYWsiPcKoXFV8xlYuwK5w8b79HYkHShp\nluf/Os+3tZU02rfNknSibz/K83TNkDTGU6jgZS+X9Kik1z1/2DWSXpA02swO8mnx6wFdPV/Yw5JW\nyIMiaVtJEyRNlfSgpO6Vu1QhhBDyVbMD2wQYZWZbkVat/9/cDknrAOcBuwF9gf6Shvjrdc1sCzPb\nkhSzBXCnmfU3sz6k9RiPyJxnTa/nRNJU+QtJOcC2lNTXy3QCppnZNsAE4NfZhvoCvpcA+5vZtqRA\n57MLvamIAwshhMqoZgeWnztsUGZff2C8mb3jubhuBHYCXge+JukSSXuQOj6ALZSyLs8izRzsnanr\nHksLPs4C3vbV65cAz5GecUFaAPjWIm2B1NluATzk+cB+RRq1rSDiwEIIoTKq+QxshfxemdeFYrQw\nsw8k9SElyvwp8APSyvSjgSFmNkPSMGCXzGG5nF9LWD7/1xKKv//8tgl4rsDqHCVFHFgIITSdao7A\n8nOHPZ7Z9wyws6Suktr6/gmSugJtzGwMcBqwjZfvAszzW30HU39tgP399UF5bQF4CeiWa6+k9r42\nYgghhCqp5ggsP3fY5cB3AcxsnqRTSQvrCvibmd3lo69rJeU63lP982mkTu8N0q3CbELLcnwM9JY0\nlRQLNjS703OE7Q9cLGl10nX7I+k2ZAghhCpotHxgtUzSQjPr3Nj1Rj6wEEKov8gH1gxEHFjrETFe\nIVRei+/AiuT8Ggv8CehGyvuViz0bTZrZ2A/4KnCymd3h+04iTRrpAIw1s+Wm2ocQQqisFr2Yb4mc\nX6OAYz2mawRwWeaw7qRp9HsD53o9uwO9gAGkWLRtJe1U5JwRBxZCCBXQ0kdgS3N+AUi6h7Rs1PbA\n7dLS2fodMsf81ePEnpe0tm/b3T+e9a87kzq0ifknNLNRpA6SDt17xQPGEEJoIi29AysUT9YG+NCX\nkCokGyumzOdzzOzK+pw84sBCCKHptOhbiBTO+fUJMEfSAQBK+tRRz4PAj70OJK0r6StN2fAQQgil\ntegRWImcXwcDl0v6FdAeuMXLFKtnnKTNgKf8tuNC4BDgP037DkIIIRRT1RFYhfKCFcr5NQf4K/B7\nM9vczM4EMLNhuVmH/nXnzOuLfAHhHYG/mNlrjdS+EEIIDdAsR2CNlRdMUjuK5PwysysaWO0apJXz\nL6urYMSBtQ4RAxZCdTSHZ2D1zgvm2xdK+n1+Di8/9reSJgDHk5abeg/4DBgsqYeXGylphL/eSNID\nnuvrMUmb+va1JY31PGMzJG1Pmlq/kaTpks6v4HUKIYSQ0Rw6sIbkBYPSObzWMLOdzez3wKXA9V7/\njcDFBdpQLC7sYmCC5xnbhrT24SnAa2bW18xOyq8o4sBCCKEymkMH1pC8YFA6h9etmdcDgZv89V/y\nyuEzC3NxYdOBK0nBzJA6zssBzGyxmdXZI0U+sBBCqIzm8Ays3nnByqjn43qcr664sAaLOLAQQmg6\nzWEEVu+8YL6vrhxeOU8CP/TXB+eXM7MFFI8Lexj4iW9vK2k14CPqn64lhBBCI2sOHVguL9hMYC38\nlh2kvGCkSRiPkuK0ppnZXb47m8NrN+DMIvUfBxzu9R9Kmtix9BT++WDgCEkzSM+59vHtxwO7SppF\nWgS4t5m9Bzzhk0piEkcIIVRJzeYDW9kcXpIuIXWI1zbw+NHAvdm4sXyRDyyEEOov8oGVIOksYDtg\nZFOeJ+LAWr6IAQuhehrlFqKvqPFigXiu/pKe9BiqSZK6+LOkCzy2a6akY72Ob0p61rdfI6mDb5/r\ncV1P+fT0bSQ9CLwt6Wgvs4ukiR6z9bykKyS18X0rxJGZ2WnA5n47EEn7+4gKSaMlXeztfl3S/r5d\nki71+u8DYi3EEEKoosZ8BpYfz3UMaTr78R5HNRj4FBgObAhsnYvNktQRGA0M9eWa2uGTJ9ybZjYQ\neMzL7Q98g+Wfew0A/g/YEtgI2LeOOLJSVsgJBnzf3+OWwFGkqfcriDiwEEKojMbswPLjub4NzDOz\nyZBm+3ks12DgCn+Nmb1P6hjmmNnLfvx1LIv3ArjbP88CnjGzj8zsHWCRpDV83yQze93MFgM3kzqg\nUnFkpfzVzJaY2fNALifYTsDNHg/2FvBIoQMjDiyEECqjMZ+B5c8GWcDyiSJzVKBsXfFeuRxdS1g+\nX9cSlr2HQvFkperNlu9Y5Hz5bavXjJeIAwshhKbTmCOw/Hiup4F1JPUH8Odf7YBxwNH+GklrAS8C\nPSV93Y8/lGXxXuUaIGlDf/Y1lBTvVSqO7G1Jm3n575dR/0Tgh/4Mrzuwaz3bF0IIoRE1ZgeWH891\nCakjucTjqx4ijXSuBv4BzPTtB5nZIuBw0nJOs0gjq/quFv8U6XnVbGAOMLaOOLJTgHtJtwLnlVH/\nWOAV0m3My6l/BxtCCKERNUocmKSepJioLep53GjqiKUqs55dgBFmtnfe9iHAy/4sq9TxI4GFZnbB\nyrQjX8SBhRBC/UUcWDKENMoq2YE1lYgDa/kiDiyE6mmUW4hmNrfU6MvjxF6QdJWk5ySNk7RKXpmC\n8V6SXsvFe3m5kz2ua4ak3BT3XsDavm2Mx6BtD3wPOF8pd9dGko6SNDlbrkBb+0p62mPUxkpa07f3\n921PSTpfBTJJhxBCqJxKroXYC/iTmfUGPgT2K1CmZLyXpD1Jo6rtPLbsd37cnWbW37e9ABxhZk+S\npt+f5Lm7XitUrkAbrgd+7jFqs1iWZ+xa4Ghv3+JibzLiwEIIoTIq2YHNMbPp/noq0LNAmbrivQYD\n15rZJ7A0hgxgC6VMyrNIC/P2LtKGkuUkrU5KhpmboHEdsJOfu4t3irAsv9gKIg4shBAqo5LPwLKx\nVYuBVUqUKRbvVSiGDNJobYiZzZA0DNilSBvKLZevPnnJloo4sBBCaDrNIZ1KfYwDfpx7duUxZJDy\nc82T1J40ssrJz91VrBwAnnH5A0k7+qZDgQlm9gHwkaRv+PYf5h8bQgihsmpqFqKZPSCpLzBF0mfA\n34BfAKeRgpbfIN1+zHVatwBXSTqO9DytWLmsw4ArvJN8nRSfBul52VWSPgbGA/GAK4QQqqhm84E1\nFUnjSTFlU/K2dzazhf76FKC7mR1foIqlIg4shBDqL+LASpDU1hf9rY+9JJ1KumZvAMPqOiDiwKor\nYrRCaNlq7RlYnVQ8N9lcSadLehw4oFi8lzvE84HNljTAt90LTCOlhFmHIulUQgghVEaL68Bcfm6y\n//Xti8xskJndQvF4L4BOZra9H3eNb/sl8IiZ9Sct5Hu+pE75J444sBBCqIyW2oHl5yYb5K9vheLx\nXpnjbwYws4nAah4HtjtwiqTppEkcHYEe+SeOOLAQQqiMlvoMrFBuMICPV+J4AfuZ2UvlNiLiwEII\noem01BFYfm6yx7M7i8V7ZYoMBZA0CJjv5R8EjpUk37d1E7Y/hBBCHVpqB5afm+zyAmUOIz3Hmgn0\nxddbdB9IepKUkyy3XuJZQHtSHrPZ/nUIIYQqafRbiE2VW6vAeeYC/YDOrJiLbImZHZ13SM/sF74u\n4zey23x5qdlmdmqBXGL3k+LD/qfR3kQIIYQGa6nPwBrDSucSiziw4iJGK4Swsuq8hVhHXFVXL9PP\nV7DI6SPpEUmvSDrKy+wiaYKk2yS9LOlcSQdLmuT5vTbyct08V9dk/9jBt39ZKY/Ys5KuZPkFdtvK\nc40Bo4D+fsxGkh6QNNVXod/Ut39X0jNe198lrZ33nlfIJea7DvD2vpx5fhZCCKEKyn0GViyuqpit\ngL2AgcDpktbx7X2A44EtSRMnNjazAcDVwLFe5iLgQo+32s/3QYrTetzMtialXclOYS+Wa2wUcKyZ\nbQuMAC7z7Y8D3/C6bgFOzja+SC4xgHbe3hNYPm5sqYgDCyGEyij3FmJ+XNVxdZS/y8w+BT6V9Cgw\ngNSxTDazeQCSXiOtLg8pkHhXfz0Y2Nwn+0GKw+pCitPaF8DM7pP0QeZ8K+Qak9SZtFrG7Zm6Ovjn\n9YBbJXUHvgTMKeMaANyZPUehAmY2itRx0qF7r1hoMoQQmki5HVihuKgvWDaC61hGeVgxx1c2/1eu\nLW2Agd4BLuWdULEOoVCusTbAh2bWt0D5S4A/mNndknYBRhapt9h5FlPGtYs4sBBCaDrl3kIsFFc1\nF9jWt+2XV34fSR0lfZmUNHJyPdo0Djgm94WnTwGYiOfwkrQnsOaKhy5jZguAOZIO8GMkqY/vXh34\nl78+rEgV+bnEQgghNCPldmCF4qrOAC6S9BhpRJI1CbgPeBo4y8zeqkebjgP6+YSR54HcdPgzgJ0k\nTSMt6/SPMuo6GDhC0gzgOWAf3z6SdGvxMeDdIsfeApzkEz02KlImhBBCldSZD0xST1aMs2pSfltv\nhJntXWDf34CDzOzDSrXHz7vQzDrX55jIBxZCCPWnlpoPzMy+01R1S2pnZl80Vn0tOQ4s4rhCCNVW\n5y1EM5vbGKMvjyd7IRev5TFdq0j6usdizZA0LXO7bjXP0/W8pCsktfF65krqKqmTpPv8uNmScusX\nbuvxZlMlPegzDZF0lMeVzfA4s1V9+2hJf/DZkudJ6izpWo9Nmylpv8x7ONuPfzo/diyEEEJlVXot\nxELxWjf6tj6kae/zvOwA4P9IMWMb4VPoM/YA3jKzPt7BPiCpPWmG4f4e+3UNcLaXv9PM+vt5XmDZ\nGocAGwODzez/gNNIC/hu6XFvj3iZTsDTfvxE4KhGuB4hhBAaqNIdWH681obAumY2FsDMFpnZJ75/\nkpm9bmaLSfm5BuXVNQsYLOk8STv6ivGbAFsADynl7foVKeYLYAtfjWMWaXJH70xdt/t5IMWh/Sm3\nw8xy8WafkZaWyrW9Z6E3GIHMIYRQGZV+BpYfr1VqKnyxWLL0hdnLkrYFvgOcI2kcMBZ4zswGsqLR\nwBAzm6G0aO8umX3ZPGEqcG6Az23ZjJeicWARyBxCCJVR7UkcC4B/ShpiZn+V1AFo6/sGSNoQeIOU\nn2tU9kBfnup9M7tB0kJgGHAu0E3SQDN7ym8pbmxmz5Fiuub5toNZFgeWLxeHdoKfZ83MKKxeIpA5\nhBCaTnPIB3YocJzHmD0JfNW3P0XqkGaTlnoam3fclsAkv1X4S+A3ZvYZsD9pMsYMYDrpuRqkZ1vP\nAA8BL5Zoz2+ANX1iyAyWLXEVQgihGakzDmyFA4rEhUm6mrQ8U4PTj+TVNwzoZ2bHqIE5xiStQYoZ\nu8y/Xge42Mz2b4w21iXiwEIIof7KjQNrtBGYmR3ZWJ1XI1qDzMr5ZvZWpTqvEEIITauhz8DaSboO\n2Bp4GfgR8DfS6hlTJB0I/II0IeI+M/s5gKQ9gN+SnnO9a2bflLQWabr714BPgOFmNrPYiT1O7E9A\nNy9/lJm96HFZV3g9AD8hLUu1kd9mfMiPu9fMtpDUkbQkVj/SwsQ/M7NHfeT3PWBV0vT9sWZ2sqS2\nwJ+9vAHXmNmFpS5SrQYyR5ByCKEWNLQD2wQ4wsyekHQNmVGO36Y7j7TQ7wfAOElDgCeAq4CdzGyO\nd1yQ1jh81syGSNoNuB4otIJ8zijgaDN7RdJ2pBxfuwEXAxPM7Pve2XQGTgG2yK1I77c/c34KYGZb\nKiW6HCdpY9/Xl9Q5/xd4SdIlwFdIU/638LrWqN8lCyGE0Jga2oGVyg/WHxhvZu8ASLqRlMtrMTDR\nzOYAmNn7Xn4Qvpq9mT2ilHl59UInVekcX7uRRoJ4TNd8SaWm6Q8iBT3jI7g3SAHNAA97XBm+oPAG\npMWAv+ad2X0sy2WW38bhwHCAtqt1K3H6EEIIK6OhHVipGC1RWLH4qkLli80sKZXjq76KtRNWjFdr\nZ2YfeDqWb5NGbz8Afpx/YMSBhRBCZTS0A+uRi7ViWX6w7/q+Z0hpVrqSbiEeSBrpPAX8SdKGuVuI\nPgrL5fk6y1ehf9fMFmRGWEv59jmSDjCz25UKbWVmM4CHSc+9/ui3EDtROqdX7ryP+K3DHsBLwDaF\nCvv7+czMxihlkx5d10WKOLAQQmg6DZ2FWCg/GABmNg84FXgUmAFMM7O7/JbicOBOj6+61Q8Zief/\nIsV9FUswmVMsx9fxwK6+VNRUoLeZvQc84TFd5+fVcxnQ1svfCgwzs/9S3LrAeJ8QMtrfYwghhCqp\ndxxYvU9QobixIuceArxcren9EQcWQgj1V/E4sPqqUNzYEGDzJj5HCCGEKqjUCOwB0rOxsuPGCsVd\nAXeTVo7fxuvuBdxiZttKOpcUv/UFaYbgnaTV4+f7Ry6vV6EYstHAp8CmpBmHh5NuZQ4EnjGzYQ2J\nA+vQvZd1P+yPDbxy1RNxYCGEaip3BFapxXwbEjf2JnlxV2b2oaT5kvp6WpbDgdEeU/Z9YFMzs0zZ\nu0m3L+/wOh6mcAwZpJXxdyN1gvcAOwBHApMl9SUFX0ccWAghNBOVuoWYHzeWze21NG7MzL4gJbjc\nCXgdj7vyFTwWePmrgcN9RDQUuMn3LQKulrQvaXS1nLwYsunAlUD3TJF7PF3KLOBtM5tlZktIE0V6\nlmhP/nkiH1gIIVRApTqweseNeQqTPsB4UtzV1b5rDLAnsDcw1cze845vgO8bQrplmW9pDFnmY7PM\n/twMxCUsHwe2BI8DK9Ke/HaPMrN+Ztav7aoF47FDCCE0gkrdQqx33FixuCszWyTpQdLU/SNg6ehq\nVTP7m6SngVe97qVxYHXEkNUp4sBCCKF5qdQIrN5xY5SOu7qRNIrLLefUBbjX658AnOjbbwFOkvSs\nLwJcLIasHBEHFkIIzUiTz0JsCpJGAKub2WnVbkspEQcWQgj119xmITYaSWNJaU52q6tsA+r+HrC5\nmZ3b2HWHEEJoXDU5AqsV1YwDi1iuEEKtavYrcTQlST/z9Q9nSzpBUk9JL0q62rfdKGmwpCckvSJp\ngB83TNKl/nq0pMslPSrpdUk7S7pG0gse+BxCCKGKWlwHJmlbUoDzdsA3gKNIQcpfBy4CtiKtuHEQ\nKR5tBGkVkEJywc0nkoKbLwR6A1t6cHOh80ccWAghVECL68BIndJYM/vYzBaSlpTaEZiTF5z8cCZw\nuWeRuuoKbl5BxIGFEEJl1NwkjjIUS1SZH5ycDVwudh1KBjfX1ZCIAwshhKbTEkdgE4EhklaV1Im0\nRuJjVW5TCCGERtbiRmBmNs0nWUzyTVeTVvgIIYTQgsQ0+jIUS8pZlwhkDiGE+mvV0+hDCCG0fC2y\nAysn7ktSJ4/rmuxrJe7jx/aWNEnSdEkzPWkmQFtJV0l6TtI4SatU8S2GEEKr1yI7MFdX3NcvgUfM\nrD+wK3C+T/o4GrjIzPqSsi//0+vrBfzJzHoDH7Isw/NysnFg77zzTpO9uRBCaO1a3CSOjDlmNgtA\n0tK4L0m5uK/1gO/5wsAAHYEewFPALyWtB9zp2Ztz9U33slMpEQcGjIL0DKwp3lgIIYSW3YHVFfe1\nGNjPzF7KO+4FSc8AewEPSjqSlI05W99iIG4hhhBCFbXkW4h1eRA41hNbImlr//w14HUzuxi4m3QL\nMoQQQjPTmjuws4D2wExJs/1rgKHAbE9cuSlwfZXaF0IIoYSIA2tCEQcWQgj1F3FgIYQQWrSa6MAk\nneZxXQ9JulnSCEkbSXpA0lRJj0na1MuOlnSxpCc9j9f+vl2Szve4sFmShvr2XSRNkHSbpJclnSvp\nYI8FmyVpIy/XTdIYjxubLGmH6l2REEIIzX4WoqR+pJirrUntnUaaxj4KONqnuW8HXEbK3QXQnRTz\ntSlpIsYdwL5AX6AP0BWYLGmil+8DbAa8T5pxeLWZDZB0PHAscAIppuxCM3tcUg/SJJDNCrR3ODAc\noEePHo14JUIIIWQ1+w6M1BHdZWafAki6hxSztT1wu08iBOiQOeavnrfreUlrZ+q52cwWA29LmgD0\nBxYAk81sntf/GjDOj5lFCnIGGAxsnjnfapK6mNlH2cZGHFgIIVRGLXRghfJ7tQE+9NUyCsnGbCnv\nc13li+UKawMMzHWkIYQQqqsWnoE9DnxXUkdJnUkBxp8AcyQdAEufb/Wpo56JwFBJbSV1A3ZiWcqV\ncowDjsl9IalY5xlCCKECmv0IzMwmS7obmAG8AUwB5gMHA5dL+hUpnusWL1PMWGCglzHgZDP7d27y\nRxmOA/4kaSbpuk0krZtY1NSpUxdKyl/po7XpCrxb7UY0A3Ed4hrkxHWo+xpsUE4lNREHJqmzmS2U\ntCqp4xhuZtOq3a66SJpSTixDSxbXIInrENcgJ65D412DZj8Cc6MkbU6avHFdLXReIYQQmlZNdGBm\ndlC12xBCCKF5qYVJHLVsVLUb0AzENUjiOsQ1yInr0EjXoCaegYUQQgj5YgQWQgihJkUHFkIIoSZF\nB9YEJO0h6SVJr0o6pdrtqQZJ10j6j+daa5UkrS/pUUkvSHrO19ZsdXwRgkmSZvh1OKPabaoWX0jh\nWUn3Vrst1SJpri+UPl3SSuWbimdgjUxSW+Bl4FvAP4HJwIFm9nxVG1ZhknYCFgLXm9kW1W5PNUjq\nDnQ3s2mSupAWoR7SCn8WBHTyWM72pNV1jjezp6vctIqT9DOgH7Came1d7fZUg6S5QD8zW+lg7hiB\nNb4BwKtm9rqZfUZaIWSfKrep4sxsIml1/1bLzOblYhZ90ecXgHWr26rKs2Shf9neP1rdf86S1iMt\nhXd1tdvSUkQH1vjWBd7MfP1PWuEfrbA8ST1JKYGeqW5LqsNvnU0H/gM8ZGat8Tr8ETiZtEh4a2bA\nOM/lOHxlKooOrPEVWvW+1f23GZbxRajHACeY2YJqt6cazGyxZ49YDxggqVXdVpa0N/AfM5ta7bY0\nAzuY2TbAnsBP/XFDg0QH1vj+Cayf+Xo94K0qtSVUmT/zGQPcaGZ3Vrs91WZmHwLjgT2q3JRK2wH4\nnj//uQXYTdIN1W1SdZjZW/75P6RF1gc0tK7owBrfZKCXpA0lfQn4ISkrdGhlfPLCn4EXzOwP1W5P\ntUjqJmkNf70KKTnsi9VtVWWZ2almtp6Z9ST9TXjEzA6pcrMqTlInn9CEpE7A7kCDZypHB9bIzOwL\nUt6wB0kP7W8zs+eq26rKk3Qz8BSwiaR/Sjqi2m2qgh2AQ0n/bU/3j+9Uu1FV0B141FMRTSY9A2u1\n08hbubWBxyXNIOVjvM/MHmhoZTGNPoQQQk2KEVgIIYSaFB1YCCGEmhQdWAghhJoUHVgIIYSaFB1Y\nCCGEmhQdWAghhJoUHVgIIYSa9P89OHA5yZT3KQAAAABJRU5ErkJggg==\n",
      "text/plain": [
       "<matplotlib.figure.Figure at 0x7fbe85968c88>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "df = pd.DataFrame(data=data, index=index, dtype=float)  # note the cast into floats\n",
    "df.sort_values(0, ascending=True, inplace=True)\n",
    "\n",
    "_ = log10(df).plot(kind='barh', legend=False, logx=False, \n",
    "            title=\"Number of H5N1 entries in each database\")"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {
    "collapsed": true
   },
   "outputs": [],
   "source": []
  }
 ],
 "metadata": {
  "kernelspec": {
   "display_name": "Python 3",
   "language": "python",
   "name": "python3"
  },
  "language_info": {
   "codemirror_mode": {
    "name": "ipython",
    "version": 3
   },
   "file_extension": ".py",
   "mimetype": "text/x-python",
   "name": "python",
   "nbconvert_exporter": "python",
   "pygments_lexer": "ipython3",
   "version": "3.5.3"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 2
}
